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ABSTRACT 


We  propose  , a  class  of  bounded  influence  robust  regression  estimators  with 
conditionally  unbiased  estimating  functions  given  the  design.  Optimal  estimators 
are  found  within  this  class.  Applications  are  made  to  generalized  linear  models. 
An  example  applying  logistic  regression  to  food  stamp  data  is  discussed.  ' 
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1.  INTRODUCTION 


In  this  paper  we  study  robust  estimation  in  a  model  with  explanatory 
variables  X  :  (X^.Y^)  are  assumed  to  be  independent  and  identically  distributed 

with  joint  distribution  P0(dy  |x)F(dx) .  We  consider  M-est imators  defined 
implicitely  by  the  equation 

n 

2  +(Y,.X..e  )  =  0.  (1.1) 

i=l  1  1  n 

In  (1.1),  8  belongs  to  a  subset  fl  of  !RP  and  ^  takes  values  in  IRP.  Usually.  is 

required  to  be  Fisher-consistent,  i.e., 

JJ  *(y.x.0)  P0(dy|x)F(dx)  =  0  for  all  6.  (1.2) 

which  implies  consistency  in  the  usual  sense  under  weak  regularity  conditions.  In 
this  paper  we  require  a  stronger  form  of  (1.2): 

S  ^(y.x.8)P0(dy Jx)  =  0  for  all  x  and  all  0,  (1.3) 

which  we  will  call  conditional  Fisher-consistency. 

In  the  linear  model  with  symmetric  errors,  essentially  all  Fisher-consistent 
estimators  which  are  optimal  in  some  sense  automatically  satisfy  (1.3).  This  is 
not  the  case  with  asymmetric  errors  or  for  generalized  linear  models.  Stefanski, 
et  al.  (19S6)  have  investigated  robust  estimators  satisfying  (1.2)  (with  F  the 
empirical  distribution  function  of  the  {X.})  but  not  (1.3).  However,  conditional 
Fisher  consistency  is  an  appealing  concept  because  it  does  not  involve  the 
distribution  of  the  explanatory  variables  (X^),  which  is  independent  of  the 
parameter  of  interest.  Moreover,  these  estimators  have  the  advantages  of  being 
computationally  simpler  in  certain  cases  (Section  3)  and  less  affected  by  the 
estimation  of  nuisance  parameters  (Section  4). 

Recall  some  general  results  and  definitions  from  robust  statistics  (see 
Hampel  et.  al.,  1986).  The  influence  function  of  an  M-<_stimator  is 

icyy.x.8)  =  D~*(0)  *(y.x,0); 


(1.4) 


■4.  *mi  m  u  i  m  *.«•*  »  »,>  i 
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where 


Djo)  =  -  M  H  'Hy.x-P)  Pfl(dy|x)F(dx)L_fl. 


(1.5) 


Under  regularity  conditions.  N  (0^  -  0)  is  asymptotically  normal  with 


covariance  matrix 


where 


V0)  V8)  V9)  *  =  v(s/,)  =  V(<K0)> 

w^(0)  =  E[sK0)^(0)T]  . 


(1.6) 

(1.7) 


Finally  we  consider  the  self-standardized  influence 


.  .2  <^IC+> 
s(^)  =  sup  sup  _ 


V-  X  -1 

sup  sup  _  =  sup  4>  W 

y.x  xTy(+)x  y.x 


(IS) 


which  measures  the  maximal  influence  an  observation  can  have  on  a  linear 
combination  of  interest  standardized  by  the  standard  deviation  of  this  linear 

9 

combination.  Integrating  (1.8)  shows  that  s(^)  >  p.  For  other  measures  of 

influence,  see  for  example  Giltinan.  et.  al  (19S6).  The  usual  goal  of  bounded 
influence  estimation  is  to  minimize  the  asymptotic  covariance  (1.6)  subject  to  a 
bound  on  the  measure  of  influence,  in  this  case  (1.8).  As  in  Giltinan,  et.  al 
(10S6),  optimal  estimators  satisfying  this  goal  depend  on  the  measure  of 
inf luence . 


2.  l.IMFAR  REGRESSION  WITH  ASYMMETRIC  FURORS 


The  purpose  of  this  section  is  to  show  that  the  differences  between  (1.2)  and 
(11)  are  relevant  also  in  linear  regression  if  the  errors  are  not  symmetric.  We 


consider  the  model 


h  *  8o  *  \9  *  V 


(2.1) 


where  (u.)  are  independent  and  identically  distributed  with  common  density  g 
having  no  point  of  symmetry.  The  usual  design  matrix  is  assumed  to  be  of  full 
rank.  The  vector  0  is  of  length  (p-1).  The  must  important  proposals  for 


.••'A/.  \  aO.A'S  ‘.''■.'aV 


-  3  - 


M-estimators  in  linear  regression  are  the  Mallows  and  Schweppe  forms: 

*u(y.x.fl)  =  W(x)  ^(y-eo-xT0)(l  XT)T:  (2.2) 

V'gfy.x.e)  =  w(x)  *>((y-0Q-xT0)  /  w(x))(l  XT)T.  (2.3) 

Here  w(x)  is  a  scalar  weight  function  and  ^  is  a  scalar  function.  Without 
additional  assumptions  on  g,  0q  is  not  identifiable,  but  the  following  result 
gives  conditions  for  consistency  of  0.  For  related  results,  see  Carroll  (1979). 

THEOREM  2.1  : 

i)  For  the  Mallows-form  (2.2),  there  is  a  constant  r  such  that 

p(y .x,0q+t,0j . 0p)  satisfies  (1.3). 

ii)  If  the  x^  are  symmetrically  distributed  with  center  c  €  |RP  and  if  w(X-c)  = 

W(-X+c),  then  there  is  a  constant  t  such  that  for  the  Schweppe  form  (2.3) 
p(y  ,x,  0o+t,0  . 0  )  satisfies  (1.2),  but  in  general  not  (1.3). 

PROOF:  i)  is  obvious  by  defining  t  as  the  solution  of  f  p( u-T)g(u)du  =  0.  In 
case  ii)  define  r  as  the  solution  of 

ff  w(x)  *>((u-t)/w(x))  g(u)  F(dx)du  =  0. 

This  is  just  (12)  for  the  first  component  of  p.  For  the  other  components  we 
observe  that 

JT  x j  w(x)  *((u-t)/w(x))  g(u)  F(dx)du  = 

//(Xj-CjJwtx)  o((u-t)/w(x))  g(u)  F(dx)du  =  0 
by  symmetry,  so  (1.2)  holds.  Equation  (1.3)  does  not  hold  in  general,  because  the 
solution  t  of 

/  «*>{(u-t)/w (x))  g(u )du  =  0 

depends  in  most  cases  on  x.  □ 

Thus,  the  Mallows  form  has  the  advantage  of  being  a  conditionally  unbiased 
estimating  equation  even  when  the  (Xj)  or  {u^}  are  not  symmetrically  distributed. 


3.  GENERALIZED  LINEAR  MODELS 


We  consider  the  canonical  form  of  generalized  linear  models,  where 

p0(dy|x)  =  exp{y  xT0  -  G(xT0)  -  s(y)}  p(dy).  (3.1) 

If  g  is  the  derivative  of  G,  the  likelihood  score  function  is 

«(y.*.0)  =  {y-g(xT0)}x.  (3.2) 

Note  that  (3.2)  satisfies  (1.3),  so  that  the  score  is  conditionally  unbiased. 

2 

Because  ?  is  proportional  to  x.  the  influence  is  unbounded,  i.e.,  s  (?)  =  00 . 

We  are  looking  here  for  M-estimators  satisfying  (1.3)  and  s (?)  <  b  which 
minimize  V(^)  in  some  sense.  In  analogy  with  a  general  principle  for  constructing 
optimal  robust  estimators  (Hampel,  et.  al  (1986).  Section  4.3).  we  consider  the 
following  ^  function: 

1 

*cond(y'x'0,B)  =  d(y.x.9.B)  wb(id(y.x.0,B)|(xTB'1x)2)  x.  (3.3) 

where 

1 

d(y ,x, 0.B)  =  y-g(xT0)  -  c(xT0,  b/(xT  B_1x)2) 

and  wb(a)  = 

where  is  the  Huber  function  H^fa)  =  min(a.b)  (a  >  0). 

We  work  within  the  context  of  the  Schweppe-form,  although  related  results  are 

obtainable  for  the  Mallows-form  as  in  Stefanski,  et  al  (1986).  The  major  change 

is  that  *b  in  (3.3)  factors  into  two  parts.  The  first  depends  only  on  x  and  is  of 
T  -1  1/2 

the  form  w^((x  B  x)  ).  The  other  depends  only  on 

d(y,x,0)  =  y  -  g(x^0)  -  c(x^0,b/(x^B  ^x)^/2) 
and  is  of  the  form  Wg( |d(y ,x. 9) | ) . 

The  function  c  and  the  matrix  B  in  (3.3)  will  be  chosen  so  that  the  side 
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conditions  (1.3)  and  s(^concj)  =  b  are  satisfied  By  the  definition  of 


(1.3)  holds  if  and  only  if  for  all  0  and  all  a  t- . 


J( y-g(P)~c(P.a))  wa(  |y-g(p)-c(P,a)  l)exp(y  0~G(p)-S(y))  pfd? 


(3-4) 


First  we  discuss  the  existence  of  a  solurion  to  (3  4). 


LEMMA  3. 1  For  any  a  >  0  and  0.  them  is  a  solution  c  =■  c(0.a to  (3.4). 


PROOF:  For  fixed  y.  0.  a,  the  funr'ian 


c  - »  (y  -  g(0)  -  c)  wa(iy  -  g (.0)  -  c|) 


is  continuous,  bounded  and  monotone  nonincreasing  with  limits  ±  a.  Hence  the 


existence  follows  from  dominated  convergence  and  the  intermediate  value  theorem. 


A  practical  advantage  here  is  that  often  the  function  c  can  be  calculated  in 


closed  or  almost  close  form.  This  is  particularly  important  compared  to  (1.2) 


where  c  is  a  vector  (depending  only  on  0)  whose  computation  is  quite  difficult. 


see  Stefanski,  et.  al  (1986),  Section  2.4.  Here  are  two  examples  where  c(0,a)  can 


be  calculated  explicitly. 


tlstic  Regression  Here  p  puts  equal  mass  at  0  and  1,  S(y)=0,  and 


G(0)  =  £n{l  +  exp(0)}.  Write  p  =  exp(0)/(l  +  exp(0))  and  q  =  1-p.  It  is  easily 


checked  that 


fap/q  -  p  if  0  <  0  and  a  <  q 

q  -  aq/p  if  0  >  0  and  a  <  p 


otherwi se 


satisfies  (3.4). 


L 2  :  Negative  Exponential  Regression  Here  p  is  Lebesgue-measure  on 


[0,®),  G (0)  =  -Bn(-p),  s(y)  =  0  and  0  <  0.  Two  cases  occur. 


If  the  bound  a  is  large,  the  Huberization  in  Is  one-sided  (for  large 


-  - 


"  ^  *  ’  v*  V 


1 


y's  only),  for  small  a's  both  large  and  small  y's  will  be  Huberized.  It  can  be 

checked  by  straightforward  calculations  that  the  cutting  point  between  the  two 

20a 

cases  is  given  by  the  equation  e  K  =  1  +  /3a,  so  that  0a  ~  -  0.797.  In  the  former 
case  c(0,a)  =  -0  *  times  the  smaller  solution  of  exp(x+0a-l)  =  x  and  in  the  latter 
case  c(0,a)  =  -0  1 ( 1  +  log(0a/(exp(0a)  -  exp(-0a)))). 

Turn  now  to  the  matrix  B.  It  follows  from  the  definition  of  ^con(j  that 


s(^c0nd)  =  b  provided 


Eel>cond(y-,<'9'B)  +cond(y’’‘'(,’B)  J  =  B' 


(3.5) 


Equation  (3.5)  is  used  to  define  B=B(9).  B(0)  depends  also  on  the  distribution  F 

of  the  explanatory  variables.  A  necessary  condition  for  (3.5)  is  b  >  p,  but  we  do 
not  know  if  it  is  also  sufficient. 

We  have  the  following  optimality  result  for  ^con(j- 
THEOREM  3 . 1  :  Suppose  that  for  a  given  b  (3.5)  has  a  solution  B(0).  Then  '^concj 
minimizes  tr(v('^)^('^concj)  *)  among  all  ^  satisfying  (1.3)  and 

sup  ic  v(*  r1  ic  s  b2. 

(y.x) 


Theorem  3.1  is  a  corollary  of  the  following  analogy  to  Theorem  1  of 
Stefanski,  et  al.  (1986).  Note  that  Theorem  3.2  below  also  applies  to  any  kind  of 
model  with  explanatory  variables. 

THEOREM  3.2  :  Let  5(y,x,9)  be  the  likelihood  score  function.  Define  the  score 


function 

1 

*cond(y-X>0)  =  (^"•in^l.b^Uf-cjVV-c)}). 
where  c  =  c(x,0)  and  B  =  B(0)  are  assumed  to  exist  and  satisfy 

E(*cond(y’Xl9)|x)  =  ° 

E<*cond(y'X-0)  ^cond^y-X-0)T}  =  B' 

Then  (3.6)  minimizes  tr (v(^)V(^cond)  l)  among  all  4>  satisfying  (1.3)  and 


(3.6) 


-  7  - 


Sup  IC  v(+  )-‘  IC  $  b2. 

(y.x) 

With  the  exception  of  nail tiplicat ion  by  a  constant  matrix.  ^con(j  Is  unique  almost 
surely.  □ 


PROOF  OF  THEOREM  3.2  The  proof  is  almost  identical  to  tha  of  Theorem  1  in 
Stefanski,  et  al.  (1986),  once  one  notes  that  for  any  conditionally  unbiased  score 
function  ■p , 


E  c(x, 0)  p(y,x,9) 

=  E{c(x,0)  E(p(y.x,9) |x) }  =0.  O 

The  computational  complexity  of  the  conditional  unbiased  estimator  is  not 
particular  to  the  model  (3.1).  For  instance,  if  we  have  a  generalized  linear 
model  with  arbitrary  link  function  h,  we  have  to  replace  in  (3.3)  d(y,x,0.B)  by 
h-(xT0){y  -  g(h(xT0))  -  c(h(xT0).b/((xTB‘1x)1/2  |  h’(xT0)|))}, 
where  c(P,a)  is  still  defined  by  (3.4). 

In  applications,  the  distribution  F  of  the  {x^}  is  unknown.  It  is  common  to 
replace  F  by  its  empirical  distribution.  From  (3.3)  and  (3.5),  this  means  that  we 
solve 

£w h-vVV  *  °-  13  7> 

-1  ^  T  T  ~  T  1/9 

N  jfj  XiXi *  V(X1  %'  b/(xi  V  xi)  )  =  V  (3. .8) 

where 

v(P.a)  =  J( y~g(P)-c(/3,a))2  w2(y ,(3,a)exp(yP-G(P)-S(y) )  p(dy),  (3.9) 

w(y.p,a)  =  Ha(|y-g(P)  -  c(P,a)[)  /  |y-g(0)  -  c(p,a)  |  (3.10) 

In  many  applications,  improved  protection  against  outliers  through  higher 
breakdown  points  can  be  achieved  by  the  use  of  redescending  p  functions,  see 
Rousseeuw  (1984).  In  equations  (3.3),  (3.4)  and  (3.10)  the  Huber  function 
could  be  replaced  by  any  of  the  redescenders  such  as  Hampel’s  three  part  function 
p  or  the  Tukey  biweight.  The  calculation  of  c(/3,a)  is  of  the  same  complexity  as 


with  the  Huber  function.  The  breakdown  properties  of  such  estimates  remain  to  be 
studied. 


4  :  THE  EFFECT  OF  ESTIMATING  THE  MATRIX  B 


In  the  last  section  we  derived  the  estimator  defined  by  (3.7)-(3.8)  as  an 
approximation  to  the  optimal  estimator  which  uses  ^^^(y.x,  0,B(0)) .  We  may 
consider  (3.7)  and  (3.8)  as  an  M-estimator  for  both  9  and  a  nuisance  parameter  B. 
The  ^function  defining  this  M-estimator  is 

^cond(y,X,S'B)T'  x(x.9.B)T)T  •  where 

T  T  T  -1  1/9 

X(x,8,B)  =  x  x  v(x  0,  b/(x  B  x)  )  -  B. 

The  influence  function  of  this  estimator  is  (compare  (1.4)  and  (1.5)) 

IC^  x(y,x,0.B)  =  D"|x(0)  (*cond(y.x.9.B)T.  \(x,0.B)T)T.  (4.1) 

where 


D 


'P.  x 


'a/3  E0 


'a/3  E0 


'/,cond(y-X'P’B) 


/3=0 


x(x.P.B) 


/3 =0 


"3A  E0  *cond(y,x,9,A) 


A=B( 0) 


A=B( 0) 


-gA  Eg  |x(x.0.A) 


(4.2) 


By  the  definition  of  'PCQn<^  and  c(/3,a)  in  (3.3)  and  (3.4).  ^cond(y  • x*  A)  satisfies 
(1.3)  for  arbitrary  A.  Hence.  Eg['/'con^(y,x,0,A)]  =  0  for  all  A  and  the  upper 


right  block  of  D,  is  zero.  This  means  that  the  9  part  of  the  influence  function 

P,x 


for  (3.7)  and  (3.8)  is  equal  to 


<-  f/e^cond l^)'1  +cond^-’<'8-B)- 


(4.3) 


On  the  other  hand,  the  influence  function  for  the  optimal  '/'o(y.x,0)  = 


,(y,x,0,B(0))  is  also  equal  to  (4.3)  because  by  the  same  argument. 


I'LlU'U't 


K'rrrrrjr?  rr 


D*  (0)  = 

M3 


op  p=Q  3A 


I  ^-B(0) 

A=B  30 


We  have  thus  shown 


THEOREM  4. 1  •’  The  0  part  of  the  influence  function  in  the  case  that  0  and  B  are 
simultaneously  estimated  by  (3.7)  and  (3.8)  is  the  same  as  the  influence  function 
in  the  case  that  0  alone  is  estimated  using  the  optimal  ^(y.x.O)  = 

A 

^cond^,x’^’®^^  ‘  Part^cu^ar'  the  asymptotic  covariance  matrix  of  0^  is  the 

same  in  both  cases.  □ 


REMARKS  : 

i)  0n  and  B^  are  not  asymptotically  independent:  EeC^concj  X^]  =  0  by  (1.3),  but 

—  Efl[x(x,P.B)]  n  /  0  in  general. 

dp  P 

ii)  Because  in  linear  regression  with  symmetric  errors  \  does  not  depend  on  0,  an 
analogue  to  Theorem  4.1  is  obvious.  In  addition,  estimation  of  the  scale  of  the 
errors  does  not  change  the  asymptotic  covariance  either,  and  0^  is  asymptotically 
independent  of  all  nuisance  parameters. 

iii)  From  the  finite  sample  interpretation  of  the  influence  function,  (4.3)  means 

A 

the  following:  to  the  first  order  of  approximation  the  change  in  0^  caused  by 

adding  or  deleting  an  observation  at  (x,y)  is 

-"Kl  ^  1  w^’VV- 

-A  A 

i.e.  the  change  in  B^  has  approximately  no  effect  on  the  change  in  G^.  In  this 
sense  the  estimator  (3.7)-(3.8)  is  reasonably  stable. 

(iv)  For  the  Fisher-consistent  estimator  (2. 12)-(2. 13)  of  Stefanski,  et  al. 

(1986),  there  is  no  analogue  of  Theorem  4.1.  The  0  part  of  the  influence  function 

T 

is  in  general  a  linear  combination  of  'J'qt.  I  x3  and  ^aC'^'rt'^rt  !x1_B  because 


'0LrBrBI 


all  blocks  in  D  are  in  general  different  from  zero. 


5:  EXAMPLES 

We  illustrate  the  conditional  Fisher-consistent  estimators  on  two  examples  of 
logistic  regression,  the  first  considered  by  Stefanski,  et  al.  (1986). 

Example  5. 1  :  The  Food  Stamp  Data 

The  response  is  whether  or  not  one  participates  in  the  federal  food  stamp 
program.  Two  dichotomous  predictors  are  available,  tenancy  and  supplemental 
income.  An  additional  predictor  is  log(l  +  monthly  income).  There  were  150 
observations,  with  24  participating  in  the  program.  One  observation,  #5,  is  known 
to  be  highly  influential  on  the  ordinary  analysis,  with  another,  #66,  slightly 
less  influential.  We  used  the  Huber  function  with  two  values  of  b,  and  in  one 
case  contrasted  the  use  of  the  conditionally  unbiased  score  function  and  the 
biased  score  function  with  c(j3,a)  =  0.  The  latter  choice  was  used  by  Stefanski, 
et  al.  (1986)  because  of  numerical  problems  in  enforcing  unconditional  Fisher 
consistency  (1.2).  In  Table  1,  we  list  the  results  of  this  analysis.  We  also 
list  the  analysis  one  would  obtain  if  one  used  a  Hampel  4*  function  with  bend 
points  (6.14,32),  see  Hampel,  et  al.  (1986,  pages  66-67).  For  computational 
convenience,  rather  than  solving  (3.4),  in  defining  the  function  c(P.a)  we  have 
chosen  to  use  the  same  formula  as  in  Example  3.1,  which  should  not  affect  the 
results  too  severely. 

In  this  example,  the  major  difference  is  not  among  the  biased  and  conditional 
Fisher-consistent  estimates,  but  rather  among  the  choices  of  b  (maximum  likelihood 
corresponds  to  b  =  <*>).  With  decreasing  b  the  importance  of  supplemental  income  as 
well  as  the  weights  for  cases  #5  and  #66  decrease  and  the  importance  of 
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log( 1+monthly  income)  increases. 

Example  5.2  :  Skin  Vaso-constrlction  data 

These  data  have  appeared  elsewhere  in  the  context  of  robustness,  see  Pregibon 
(19S1).  The  response  is  the  occurrence  of  vaso-constriction  in  the  skin  of  the 
digits  and  is  regressed  on  the  logarithms  of  the  rate  and  volume  of  air  inspired. 
In  our  work,  we  took  Rate^  =  0.30,  see  Pregibon  (1981). 

This  data  set  is  inherently  unstable.  As  previous  authors  have  shown,  once 
one  deletes  observations  #4  and  #18.  almost  perfect  discrimination  is  possible. 
Our  analyses  are  listed  in  Table  2.  We  used  the  Huber  weight  function  only,  and 
have  selected  various  values  of  b.  When  we  tried  to  use  the  Hampel  ^  function, 
observations  4  and  18  were  immediately  given  weight  0,  in  which  case  even 
computing  the  maximum  likelihood  estimates  is  delicate. 

A  biased  analysis  uses  a  Huber  function  wi th  b  =  6.408,  but  with  c(j3,a)  =  0. 
When  we  used  this  value  of  b  in  our  conditional  Fisher-consistent  score,  we  find 
that  observations  4  and  18  are  barely  downweighted  and  the  resulting  analysis 
looks  very  much  like  the  usual  likelihood  analysis.  As  we  move  b  to  5.5  and  5.0, 
observations  4  and  18  are  given  very  little  weight,  and  parameter  estimates  and 
standard  errors  change  dramatically.  In  this  example  too  the  value  of  the 
sensitivity  b  seems  to  be  most  important,  and  we  recommend  to  successively 
decrease  b  and  see  how  estimates,  standard  errors  and  weights  change. 


6  :  CONCLUSIONS 


Conditionally  unbiased  score  functions  are  appealing  because  their  definition 
does  not  depend  on  the  distribution  of  the  predictors  {x.}.  In  the  context  of 
robustness,  the  resulting  estimators  have  an  analogous  optimality  theory  to  that 
already  developed  for  unconditionally  unbiased  score  functions.  In  addition, 
conditionally  unbiased  score  functions  are  often  far  easier  to  define.  Although 
ignoring  the  bias  and  setting  c  =  0  turned  out  rot  to  matter  much  in  the  examples 
considered.  one  can  construct  situations  where  this  bias  is  large.  With  our 
estimator,  we  avoid  this  problem  with  little  additional  complexity. 
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TABLE  1 


This  is  a 

reanalysis 

of  the  food 

stamp  data. 

For  selected 

observations,  the 

weights  w,  in 

D 

equation 

(3.3)  are  computed. 

Huber (7) 

Huber(7) 

Huber (5 . 5) 

Hampe 1(6.14,32) 

MLE 

c(P,a)=0 

Conditional 

Cond i t i ona 1 

Condi t ional 

Unbiased 

Unbiased 

Unbiased 

Intercept 

0.93 

4.26 

4.51 

5.49 

6.00 

(1.62) 

(2.55) 

(2.54) 

(2.66) 

(2.76) 

Tenancy 

-1.85 

-1.85 

-1.78 

-1  .76 

-1.80 

(.53) 

(.54) 

(.54) 

(.51) 

(  54) 

Supplemental 

Income 

0.90 

0.75 

0.74 

0.62 

0.70 

(.50) 

(.52) 

(.51) 

(.52) 

(.52) 

Log( 1+MI ) . 
MI=Monthly 

Income 

-0.33 

-0.89 

-0.93 

-1.10 

-1  .  18 

(.27) 

(.43) 

(.43) 

( .45) 

(.47) 

Weights 

#5 

0.21 

0.  16 

0.  13 

0.0 

#66 

0.76 

0.60 

0.41 

0.54 

For  selected 


This  Is  a  reanalysis  of  the  skin  vasoconstriction  data, 
observations,  the  weights  in  equation  (3.3)  are  computed. 


MLE 


Huber{6. 41 ) 
c(/3.a)=0 


Huber (6 . 41 ) 
Conditional 
Unbiased 


Huber (5. 5) 
Condi tional 
Unbiased 


i 

■ 

w 

i 


Intercept 

-2  92 

-5.71 

-2.98 

-6.41 

(1  29) 

(2.45) 

(1.35) 

(2.84) 

1 og( volume) 

5.22 

9.  13 

5.27 

9.98 

(1.93) 

(3.73) 

(1.93) 

(4.3S) 

log( rate) 

4  63 

8.09 

4.67 

8.S5 

(1.79) 

(3.31) 

(1.S6) 

(3.S2) 

Weights 

#4 

0.38 

>.S0 

0.25 

#18 

0.44 

>80 

0.29 

■S' 
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